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,__! ' Abstract 

^o , A one-dimensional white-in-time passive scalar model is introduced. Strong and 

persistent structures are shown to be present. A perturbative expansion for the 
scaling exponents is performed around a Gaussian limit of the model. The resulting 
predictions are compared with numerical simulations. 
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6. 

It is commonly observed in turbulent flows that coherent and well localized structures 

emerge in the sea of disorganized, random regions. Well-known examples are provided 

by 2D and 3D Navier-Stokes turbulence : numerical simulations show that structures 

appear respectively in the form of quasi-circular vortices and thin, elongated vortex tubes. 

No ab initio theory for the statistical properties of the flows (and a fortiori of their 

dependence on the structures) is however available. A systematic statistical theory has 

recently been developed for white-in-time passive field models. Their crucial property is 

that simultaneous correlation functions obey closed equations of motion [|l| . Anomalous 

scaling and intermittency have in particular been considered. The general mechanism for 

the appearance of anomalous scaling is associated with zero modes of the equations for 

the correlation functions || ]3|, |J. For Kraichnan passive scalar model, anomalies first 

appear at the level of fourth-order correlations [1-3, 5]. When a large-scale gradient is 

imposed, also the third-order moment is anomalous 0. Scaling exponents have been 

calculated by using perturbative expansions around Gaussian limits 0, |3|. For magnetic 

fields a nonperturbative anomalous solution for second-order correlations has been found 

M, The zero mode responsible for anomalous scaling comes from the balance in the 

inertial range between field lines stretching and eddy-damping. This suggests that also 

for white-in-time models, structures could be present and play an important role. A crucial 



point is whether the perturbative expansions around Gaussian limits are capable, as those 
in critical phenomena, to correctly capture the global statistical effects of structures. 
Comparing the perturbative predictions with the results of numerical simulations would 
clearly help clarifying this issue. The 5-correlation in time of the advecting velocity makes 
however simulations quite intricate ||. We have therefore concentrated our attention on 
the ID white-in-time passive scalar model presented here. As in Kraichnan model, the 
velocity is Gaussian and 5-correlated in time. Its structure function scales with a positive 
exponent £ and the dynamics has a Gaussian limit for £ — > 0. 

The simulations show that the activity of the scalar is indeed concentrated in strongly 
localized peaks. Despite of the ^-correlation in time of the advecting velocity, their lifetime 
can be significant. The mechanism of formation of the peaks is associated with the 
stretching due to velocity gradients. The peaks are responsible for the observed very high 
tails of the probability distribution functions (p.d.f.). 

The 5- correlation in time ensures that correlation functions satisfy closed equations 
of motion. The behaviour of second-order correlations is obtained analytically and their 
scaling is normal. Predictions for the scaling exponents of higher-order correlations are 
derived using the perturbation scheme in £ proposed in Ref. ||. Analytic predictions (in 
the form of Pade approximants) for the fourth and sixth-order structure functions are 
compared for £ = 0.5 with numerical results. 

The equation of the model is 

d t 4> + V (v0) = «A0 + /. (1) 

Here, both the velocity v and the injection / are Gaussian and (^-correlated in time. The 
velocity correlation function is 

(v{x, t) v(x', £')> = S{t - if) [D -S(\x- x'\)} . (2) 

The structure function S scales as S(x — x') = D \x — a/p (and < £ < 2) in the range 
of scales between the ultraviolet and the infrared cutoffs Ajjy and Am (the smallest and 
the largest scales in our problem). The large-scale injection satisfies (f(x,t)f(x',t')) = 
S(t — t') Fl(x — x') with its Fourier transform Fl concentrated around wave numbers 
0(1/ L). It is assumed that Fl vanishes at k = 0, i.e. f_~ Fl(x)cIx = 0. This condition 
is associated with the fact that is a "gradient field". The equation for the gradients 
in Kraichnan model indeed coincides with ([!]) when vector components are omitted. The 
gradient nature of appears even more clearly from the solution for the second-order 
correlations presented later. The homogeneous part of the equation for the integrated 
field 9(x,t) = f x 4>(y,t) dy (scaling with positive exponents) is 

d t e + vV6 = KA9. (3) 
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Figure 1: The gradient field <ft v s the spatial coordinate x for £ = 0.5. 



Note that for (^|) the maximum principle holds. The stretching process due to the gradients 
of v operates on the gradient field 0. To fully resolve the structures we have then preferred 
to deal directly with ([]]). This equation has been integrated by using a pseudo-spectral 
code with periodic boundary conditions. At each time step a new realization of both v 
and / is generated. The multiplicative term in (]T|) is numerically treated explicitly (a 
la Ito). The term differing in the formulations Ito and Stratonovich (the latter being 
relevant for (|l|)) is taken into account f|. The anomalous exponents are not expected to 
depend on the type of dissipation. For the measurement of the scaling exponents we have 
then used a A; 8 hyperdissipation, but we have also checked that structures are present for 
normal dissipation. The resolution is N = 2 14 . The injection is concentrated on the first 
mode. The velocity v is simply generated in Fourier space. 

A typical plot of the gradient field <ft is shown in Fig. 1. It is evident the presence of huge 
peaks of activity. Their mechanism of formation (and persistence) can be immediately 
grasped from the original equation ([]]). Neglecting diffusion and forcing, the extrema <f) 
of the field (f> obey d t <j) = — (Vt> ) <f). Since v is (^-correlated in time, the logarithm of <f) 
evolves then as a Brownian motion. On the other hand, it is known that random walkers 
have a tendency not to change sign. Specifically, let us start a walk from the origin. 
The cumulative distribution for the fraction of the walk spent on the negative side obeys 
the arc sine law [jXO | : keeping the same sign for the whole walk and equipartition of the 
time between positive and negative values are the most and the least probable events, 
respectively. In our case, this implies that there will be very long periods of time when 
the stretching mechanism operates. Peaks such as those in Fig. 1 can then be maintained 
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Figure 2: The probability distribution function (p.d.f.) of the gradient field 0, normalized 



to its r.m.s. value J {(J) 2 ). 



for significant times, despite the ^-correlation of the velocity. This is indeed observed in 
our numerical simulations both for normal and hyper- dissipation. Strong peaks persist 
for long enough to make them clearly identifiable and meaningful the use of the word 
"structures" to denote them. The structures are evidently the cause of the high tails of 
the p.d.f. in Fig. 2. 

Let us now derive the predictions for the scaling exponents. The ^-correlation in time 
ensures that correlation functions obey closed equation of motion. By using Gaussian 
integration by parts, it is for example easy to derive the equation for the second-order 
correlation C 2 (x, t) = ((f>(x, t) 0(0, t)) : 



d t C 2 = ^[(2K + S)C 2 ] + F L . 



(4) 



The solution is even and satisfies / C 2 — 0. In the inertial range of scales it decays with 
the negative exponent — £. The scaling of the second-order structure function of 6 is 
obtained by simple integration. For £ < 1, its inertial scaling is 2 — £, as in Kraichnan 
model. 

Let us now consider the scaling of higher-order correlation functions. Closed equations 
of motion can be written both for the correlations C 2n = (<f>(xi) . . . (f)(x 2n )) and those 
of the integrated field 9. Anomalous scaling is associated with the zero modes of the 
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Figure 3: The measured second-order structure function for £ 
the slope 2 — £ = 1.5. 
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homogeneous equations. The relevant equations for the zero mode Z 2n of the 2n-th order 
correlator are 
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(5) 



Fol- 



where Xjk = Xj — x k and Vj stands for the spatial derivative with respect to Xj 
lowing Refs. [[| and ||11|| , we look for a scale-invariant solution of this equation as an 
expansion in powers of £. Specifically, the structure functions are expanded as S(xjk i 

D(l+£log\x jk \) + 0(e) and Z 2 , 
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+ £, Z 2r ~ + C(^ 2 )- The equations (for any n) 
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at the lowest order in £ are trivially satisfied by taking Z 2n 
order in £, the solution is 



const. 
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The scaling exponent (at order £) of the zero mode can now be obtained operating on 
Z 2n with the Euler operator £ = X^=i ( x j^j) ■ It is easily checked that the equality 
SZ 2n = —£n(2n — 1) Z 2n holds at the first order in £. The amplitude of the correction 
to normal scaling — n£ is thus 2n(n — 1)£. This same value is obtained considering the 
equation for the correlations of the integrated field 9. Note that the anomalous factor 
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Figure 4: The measured fourth (left) and sixth (right) order structure functions for £ 
0.5. Solid lines have the slopes predicted by (^), i.e. £4 ~ 2.29 and ( 6 — 2.67. 



2n(2n — l)/2 is simply related to symmetry considerations, i.e. it is the number of couple 



12 



of points in the correlation function, as conjectured in Ref. 

For the comparison with numerical simulations, it is convenient to consider the scaling 
of the structure functions (\9(x,t) — 9(0,t)\ p ) ~ \xy p . In particular, the prediction for the 
fourth and the sixth order are 



C 4 = 4-6£ + 0(£ 2 ) 



1 + 1.5£ 



and C 6 = 6-15£ + 0(£ 2 
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l + 2.5£' 



(7) 



where we have used the classical procedure of Pade approximants |13| . For the comparison 
we have considered £ = 0.5. The measured second-order structure function is presented 
in Fig. 3. The fourth and sixth order structure functions are shown in Fig. 4. It is worth 
to note that also low order moments have anomalous scaling. The measured exponents 
for the moments of order 1/2 and 1/8 are 0.585 and 0.163. The corresponding normal 
values are 0.375 and ~ 0.094. It is an interesting issue the relation of these anomalies 
with those of correlation functions. Note also that £ = 1 separates the cases of negatively 
(£ < 1) and positively (£ > 1) correlated increments of the velocity field. We are currently 
investigating the influence of positive correlations on structures and anomalous scaling. 

In conclusion, we have introduced and investigated a ID white-in-time passive scalar 
model. Numerical simulations show that the typical configuration of the field is strongly 
structured. The possible role of structures as "atoms" for intermittency has recently 
been addressed in Ref. [I I] for shell models (for white-in-time shell models see Ref. [13]). 
The conclusion of Ref. [13] was that the single structures are not quite elementary, i.e. 
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interactions between them and with the background play a key role for scaling properties. 
This is also the case for our model. A naif argument based on the scaling of a single 
structure would indeed suggest the same asymptotic scaling as in Burgers equation, i.e. 
a constant unit value. Already the second order structure function scales however with 
an exponent larger than one for £ < 1. In the zero mode formalism used here, structures 
and their correlations do not appear explicitly, but only via their global statistical effects. 
Our results indicate that these effects are correctly taken into account by perturbative 
expansions around Gaussian limits. The agreement points in the direction of the scale- 
invariant zero mode mechanism for the intermittency of the model. 

Acknowledgments. We are grateful to E. Balkovski, B. Dubrulle, G. Falkovich, U. Frisch, 
I. Kolokolov and V. Yakhot for helpful discussions. We thank the "Meteo-Hydrological 
Center of Liguria Region" and the "Swiss Scientific Computing Center" where part of the 
numerical analysis was done. 

References 

[1] R.H. Kraichnan, Phys. Rev. Lett, 52, 1016, (1994). 

[2] M. Chertkov, G. Falkovich, I. Kolokolov & V. Lebedev, Phys. Rev. E, 52, 4924 
(1995). 

[3] K. Gawedzki & A. Kupiainen, Phys. Rev. Lett, 75, 3834, (1995). 

[4] B.I. Shraiman & E.D. Siggia, C.R. Acad. Set., 321, Serie II, 279, (1995). 

[5] A.L. Fairhall, O. Gat, V.S. L'vov & I. Procaccia, Phys. Rev. E, 53, 3518, (1996). 

[6] A. Pumir, Europhys. Lett, 34, 25, (1996). 

[7] M. Vergassola, Phys. Rev. E, 53, R3021, (1996). 

[8] R.H. Kraichnan, V. Yakhot & S. Chen, Phys. Rev. Lett, 75, 240, (1995). 

[9] P.E. Kloeden & E. Platen, Numerical Solution of Stochastic Differential Equations, 
Springer (1992). 

[10] W. Feller, An Introduction to Probability Theory and its Applications, Vol. 1, J. Wiley 
& Sons, (1950). 

[11] D. Bernard, K. Gawedzki & A. Kupiainen, Phys. Rev. E, 54, 2564, (1996). 



[12] M. Chertkov & G. Falkovich, Phys. Rev. Lett, 76, 2706, (1996). 

[13] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, Oxford Univ. Press, 
(1989). 

[14] J.L. Gilson & T. Dombre, "Towards a two-fluid picture of intermittency in the GOY 
model of turbulence" , preprint, (1997). 

[15] R. Benzi, L. Biferale & A. Wirth, "Analytic calculation of anomalous scaling in 
random shell models for passive scalar", preprint, (1997). 



